X 



Analysis of the B —> K%(— > Kiv)l + l decay 

S. Rai ChoudhurjQ 

Centre for Theoretical Physics, Jamia Millia, New Delhi 110 025, India 

A. S. CornelQ 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

G. C. Joshi and B. H. J. McKeUaJ! 

School of Physics, University of Melbourne, Victoria 3010, Australia 
(Dated: February 2, 2008) 



In this paper we study the angular distribution of the rare B-decay B — > iff (— » Kn)l + l , which 
is expected to be observed soon. We use the standard effective Hamiltonian approach, and use the 



I. INTRODUCTION 
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form factors that have already been estimated for the corresponding radiative decay B — > K 2 f- 
£S) ' The additional form factors that come into play for the dileptonic channel are estimated using the 

large energy effective theory (LEET), which enables one to relate the additional form factors to 
the form factors for the radiative mode. Our results provide, just like in the case of the if* (892) 
resonance, an opportunity for a straightforward comparison of the basic theory with experimental 
, results, which may be expected in the near future for this channel. 
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The flavour changing neutral current (FCNC) decays, like b — > sl + l~ or b — > dl + l~ , have attracted a lot of attention 
during the last decade Q, 0, this being due to the fact that these processes are forbidden at the tree level in 
the standard model (SM) but are induced by loop corrections. They are therefore very sensitive to details of the 
Oh, model and hence form a very appropriate process for the study of physics of the SM and beyond. Using the standard 
theoretical framework for the study of these processes, this framework being based on the short distance expansion 
of the relevant quark operators 1], the values of the Wilson Coefficients that enter the expansion have now been 
calculated to a very high degree of accuracy. As such the major uncertainty now lies in the evaluation of the matrix 
elements of the relevant quark operators between the physical states of the hadrons in the process being considered. 

For semileptonic decays, like B — > K*l + l~ , the relevant quark operators are bilinear in the quark fields, and the 
matrix elements of the quark operators reduce to form factors. The form factors involved in this case are those of the 
vector, axial vector and tensor current operators between the initial B and final K* meson. There have been many 
theoretical analyses of these form factors and a most promising approach for our purposes is the large energy effective 
theory (LEET) jj| . The LEET theory brings in an enormous economy of parameters since it yields relations between 
the various form factors involved 0. Once the form factors are estimated, the formalism permits the calculation 
of the full amplitude for the decay in all helicity states. Note that there have been calculations before of lepton 
polarization for these decays Q. For the if* (892) case, a calculation by Kim et a/.|?| pointed out the possibility of 
studying the azimuthal angular distribution of the decay products of the K* as a probe for new physics. Of late, 
the Babar group 0| has obtained experimental estimates of the radiative decays of the decay B — > if|(f 430)7. The 
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decay rates are comparable with the corresponding one for the if* (892). We can therefore hope that in analogy to 
the K*(892) case, data on decays of B — > K2 (1430)Z + i~ would also be availabe in the near future. We carry out 
in this work an analysis of azimuthal angular dependence of the decay products of the K* resonance similar to the 
analysis done by Kim et aZ.Q. The radiative decay for this resonance has already been analyzed by Cheng et aL0 and 
there is reasonable agreement between theory and experiment for this process. Some, but not all, of the form factors 
involved in the corresponding dileptonic mode are related to the radiative mode. However, the LEET relates the form 
factors of various quark bilinears, and these relations suffice to tackle the dileptonic mode without introducing any 
new parameters. The angular distributions involve the full density matrix rather than simply the diagonal ones for 
a decay rate, and thus would prove to be an even better test of the underlying theory than the decay rates of the 
radiative mode alone. 



II. THE EQUATIONS 



The full process, together with a representation of the kinematics we have used, is shown pictorially in figure 1. 
This figure also shows the labelling of the momenta of the various particles involved in this process. Following now 
the standard effective Hamiltonian approach, we first note that the B — > K*l + l~ process is described at the quark 
level by the process b — * sl + l~. In the limit where we neglect the mass of the strange quark, this is described by the 
effective Hamiltonian: 
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eff 
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(2.1) 



where q = (p + + As such, for the process under consideration, we need an estimation of the matrix element 
for the vector, axial vector, tensor and the pseudotensor bilinears between the \B) and \K%) states. These can be 
expressed in terms of invariant form factors as follows: 



i(K2(pK')\sa^q u j 5 b\B(p B )) = —{pb+Pk-Y ^{pk- + Ps)-q - (pk- + PB)^* (yp ( f 

r " 2 
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where q 2 = (p B - pa') 2 , q-(j>B + Pk*) = m% - m 2 K *, e* ap -PK* = °j and: 



i(K* 2 (p K »)\sa^q»b\B( PB )) = e *- Pb*P b P P k* , (2-3) 



(K 2 (p K ,)\sj p j 5 b\B(p B )} = e* Q/3 Ai (g afi p B p + gp^PBa) + A 2 (p Ba f3 P ) + A 3 (p Ba p Bf3 p K » M ) 



and finally: 



£ at Lv P P*P B PBi3 + a <-> 



(2.4) 



(2.5) 



Note that the form factors E/i, and C/3 have been estimated in the covariant light cone approach by Cheng and 
Chua 9], where their results are consistent with the data availabe for the radiative decay of the B — > K 2 j process jg. 
For the estimation of the form factors A%, A 2 , A3 and V we follow the LEET effective theory analysis given by Charles 
et oZ.|4j. As such, we obtain the relations: 
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where we have taken the limit of the heavy quark mass going to infinity and E = pb-Pk* /tub- Note that with this 
approach we have introduced no extra hadronic form factors beyond what is required for the radiative mode. Thus, 
once we are able to describe the radiative mode we have in effect a check on the model from the dileptonic mode. 

The final piece required for the study of the full decay process, B — > K 2 (-^ Kir)l + l~, involves the decay of the 
K 2 — > Kit. Note that the data for this process already exists, where the vertex is described in terms of the coupling 
constant g by: 



{K { p K )n(p^)\K 2 (p K ,)) = -^—tafiPKPK 



(2.7) 



The known width of the K 2 , T, and the branching ratio into this channel can then be used to determine the constant 
g. The amplitude for the process, as given in figure 1, can be written as: 



M = Mc +Mn , 
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Where we shall neglect the masses of the lepton, the K and the pion in comparison with the mass of the J3-meson. 
Note that this is consistent with our determination of the relations between the form factors, as we have already taken 
the heavy quark limit. 

Using the expression for the amplitude M above, we can write: 

1 2 



where Tk* 



\M\ Z = 

2 \5/2 
g A ' 

960m9-,7r' 



aGF^CKAig 
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+2 {2p_.p KP+ .p K ) \B\ 2 + 2 (2p_.p K .p + .p K , - m 2 K ,p_.p+) \C\ 2 



+Alm[DA* + CA*] 



P+-Pk(p--Pk-) 2 - P-.pk{p+-Pk-) 2 - p-.p+pK-PK- (p- -p+).pK- 



+ (P+ ~P-)-Pk {-m 2 K ,p_.p + + p_.p K ,p + .p K ,) 



\Ai 



+4Im[AB*] p + .p K *(p-.p K ) 2 - p-.p K *(p + .p K ) 2 - p-.p + p K .p K *(p- -p + ).p K 

+ {p+-P-)-Pk* (+P-.PkP+-Pk) 

+A\A\ 2 - 2p-.p K p-.p K ,p + .p K p + .p K , + p + .p-p K .p K ,(p- -p + ).p K (p_ -p + ).p K . 

+ (P-.Pk*) 2 (p+-Pk) 2 + (p+.Pk<-) 2 (p~-Pk) 2 + m 2 K ,p-.p K p + .p K p + .p- 



1}P--P+ ) ( m K*(P--PK) 2 +m 2 K ,{p + .p K f) 



4e af ^p- a p+pp Kf .PK*, {-Im[7VG* + NH*} + q.p K Re[QN*} + q.p K ,Re[QG* + QH*]} 
+4 {p- -Pk*P+-Pk + P--PKP+-PK* - P--P+pk-Pk*)R-c[NG* + NH*} 
+4 (-m 2 K »p_.p + + 2p_.p K .p + .p K .) [r C [GH*] + ±\H\ 2 ^j 
+2 (2p_.p K p+.p K ) \N\ 2 + 2 (2p_.p K ,p + .p K , - m 2 K ,p_.p+) \G\ 2 



Alm[HQ* + GQ* 



p + .p K {P--PK') 2 - P--Pk{p+-PK') 2 - P-.p+PK-PK'(P- -p+).pK* 



+ (P+ ~P-)-Pk {-m 2 K ,p-.p + + p-.p K ,p + .p K A 



-4Im[Q7V*] 



P+.PK'{p--Pk) 2 - P--Pk'(p+-Pk) 2 - P-.p+PK-PK'iP- ~P+)-Pk 



+ {P+ ~P-)-Pk> (P--PkP+-Pk) 



(2.13) 



(2.14) 



5 



+4|Q| 5 
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Note that in the above we have used the following constants: 



\m B q 2 

2k l , , A (n ( m * U - 



C = [ :-^fc 3 - 2fc!(fc 4 - hk 3 ) )[C 7L [ -H-L (m 2 B -m K ,)- A X (C 9 - do) 

■ylu u u \ _l 2 1,21, \ f 1 a in n \ n mbU ' 1 2U 3 m b n (, , m 2 K , - p b -PK- 

2{k A - k x k 3 ) + -k^ --A 3 (C 9 - do) - C 7L ^ T C 7L 1 H * 2 — 

3 J \ 2 vtibI rriBq \ m B ~ m K* 

~2k 2 \ ( fm b U 2 \ , 2 9 



G = ( ^j^fc3 - 2fc!(fc 4 - kik 3 )j \C 7L (m£ - m£.) ~ Md + do) 

2(fc 4 - fc^) + -fc£fc 5 --A 3 (G 9 + 610) - C 7L 5 -C 7L 1 H ^ 2 — 

3 J \ 2 vriBl msq \ m B — m K , 



D = \k 2 (c 7l . ) ) - -1 1 (d - f mi 1 

3 V \m B q z J 

■ylu 7 j \2 1 2 u 2u\{ 1 a in n \ n mbU ' 2 ZU 3 m b ( m%* - k, 
2{k i - kik 3 ) + -k{k 5 I I -~A 2 (d - do) - dz, _ 2 - — ^ C 7L 1 

// = -hi^C 7L \^\m%-m 2 K ,)-A 1 {C Q + C 10 ) 

2(fc 4 - kika) + -k x k 5 --A 2 {C 9 + do) - C 7L ^ tC 7 l — 5 o- 

3 / V 2 m B 9 2 \m 2 B - m 2 K 



3 J \ 2 rriBq msq \m B — m K , 

2 , 2 ( m bU2 



(2.15) 
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where k x = -4—Pk-Pk*, h = p K *-Pb, k 4 = p K -Pb and k 5 = m? - — 

"'A'* |— 1 '"a-* 

If we now use the kinematics as prescribed in Kim et al.\$, that is, where we set p = \Jp 2 K * , I = -J (j>+ + P-) 2 and 
A = \{m 2 B — p 2 — I 2 ) 2 — p 2 l 2 . Furthermore, we shall introduce various angles, namely 9k as the polar angle of the 
K momentum in the rest frame of the K* meson with respect to the helicity axis, i.e. the outgoing direction of K*. 
Similarly 9 + as the polar angle of the positron in the dilepton CM frame with respect to the K* momentum, and 
finally cf> as the azimuthal angle between these planes, that is, the K* — > Ktt and B — » K*l + l~ planes. In this case 
our constants fc, become: 
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k\ = ta — 2 — ' fc 3 = [\/\+p 2 l 2 +p , k 4 = — + — 

2m z K , tub V 12 2m b 
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and 



In which case our previously defined constants will have the angular structure: 



A = A (1) + A {2 ) cos 6< a- 
B = B (l) + B {2 ) cos9 K 
C = C(x) + C( 2 ) cos^a- 
D = L> (1 ) + D( 2 ) cos Ok + £>(3) cos 2 9 K 



Q = Q(l) + Q(2) C0S6>A 

N = N {1) + N {2 ) cos6 K 
G = G(i) + G( 2 ) cos 0a 

H = + H {2 ) COS 6<A + -ff(3) cos 2 6»A, 



where we have defined these additional constants after the various integrations in appendix B. 
The decay width for the full process can now be expressed as: 
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Note that as we have done the p 2 integration in the narrow resonance limit of K* 
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the double differential decay rate can be expressed as 
dT !->n->T7, | v tb v t 
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(2.24) 



(2.25) 



(2.26) 



Note that from now on I 2 ~ s, the dilepton frame CM energy. Again, the constants A^ etc, after the p 2 integration, 
have been defined in appendix B. 



III. NUMERICAL RESULTS AND DISCUSSSION 

The parameters in the last equation for the differential decay rate are known and summarized in appendix C. The 
fivefold differential decay rate can easily be integrated. The integration over the variable p 2 is trivial in the zero width 
approximation as stated before. Integration over the polar angles 9k and 9+ is easily done through Mathematica 
to yield the differential rate -j^- Figure 2 shows our results for this quantity which can be easily compared with 
experimental results as and when they become available. More important from the theoretical point of view is the 
dependence of this quantity on the angle 4>. This depends crucially on the nature of the Hamiltonian and also on the 
spin of the K* resonance which decays into hadrons. Following Kim et aZ.|2j we define the normalized differential 
rate: 



K<M = ^- (3-D 

ds 

We show in figure 3 a plot of this quantity, which once again is easily experimentally accessible once this process is 
seen. Our results as given in figures 2-5 show considerable structure and therefore would be a good test of the theory 
when comparing with experimental data. The results are similar to the parallel case for the similar processes with the 
K* (892) resonances. However, the two involve different hadronic form factors and thus involve different combinations 
of the Wilson coefficients. Simultaneous comparison of the theoretical results for both these channels would thus 
be a more complete check of the underlying theory. The numerical results that we have presented are for values 
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2 71 




of Wilson coefficents calculated on the basis of the SM. In theories involving new physics, which would change the 
Wilson coefficents or add other coefficients (possibly scalar and pseudoscalar hadronic currents), our expressions are 
still valid. Note that the numerical results would change depending on the exact nature of the new physics introduced, 
and thus would again be a good situtation to confront theoretical results with experimental data. In summary, we 
have presented in this paper theoretical predictions based on the standard effective Hamiltonian for FCNC processes 
together with the LEET predictions of the hadronic form factors for the process B — > K^l + l~ , which we hope will 
soon be experimentally observed. 
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FIG. 4: The ratio of the double differential decay and the differential decay rate, j^r/^, for <f> = 0. 



0.1 



0.09 



0.08 - 



0.07 



0.06 



ds d<p ds 




2 7T 



FIG. 5: 77ie ratio of the double differential decay and the differential decay rate, for s = 5GeV 2 (the solid line), ana 

s = 12GeV 2 (the dashed line). 
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APPENDIX A: THE FORM FACTORS 



The form factors we shall take from Cheng et al. 0| , where the remaining form factors can be related to these using 
the relations in Charles et aZ.0|, which lead to: 
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Therefore: 



V = iU 1 /m% , A 1 = -U 2 /2 , A 2 = , A 3 = 2U 3 /m 2 
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0T6 

1 - 2.19(s/m|) + 1.80(s/m|) 2 ' 



APPENDIX B: THE CONSTANTS AFTER p 2 INTEGRATION 

After the p 2 integration the additional constants are defined as: 
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APPENDIX C: INPUT PARAMETERS AND WILSON COEFFICIENTS 

The input parameters used in the generation of the numerical results are as follows 12]: 

m B = 5.26GeV , m K * = 1.43GeV , m b = 4.8GeV , m c = 1.4GeV , 
m s = 0.lGeV , B( J/ip(lS) -> £+£- ) = 6 x 10" 2 , 
(1S) = 3.097GeV , B(ip(2S) -> = 7.3 x 10~ 3 , 



m J/rP(lS) 

m H2S) = 3.686GeV , T H2S) = 0.277 x lQ- 3 GeV , 



V(2S) 



J/V(1S) 



0.093 x 10" 3 GeV , V tb V t * = 0.0385 , a 



12!) 



G F = 1.17 x 10~ 5 GeV" 2 . 



The Wilson coefficients used were as in Kim et al. 7], namely: 
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where the function g is taken from reference [ll| : 
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